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Note 

Based  in  part  on  the  Master's  Thesis  submitted  by  Nathan  Brilliant  to  the  Department 
of  Physics  of  the  University  of  Colorado  at  Boulder.  Thanks  to  Dana  Z.  Anderson  of 
the  University  of  Colorado  and  Judah  Levine  of  the  Joint  Institute  for  Laboratory 
Astrophysics  for  their  careful  reading  of  the  thesis  and  their  constructive  criticism. 
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Video  Microscopy  Applied  to  Optical  Fiber 
Geometry  Measurements 


Nathan  Andrew  Brilliant,  Bradley  K.  Alpert,*  and  Matt  Young 

Optoelectronics  Division 

National  Institute  of  Standards  and  Technology 

Boulder,  Colorado  80303-3328 


We  describe  a  video  microscope  used  to  measure  the  mean  diameter  and  noncircularity 
of  the  cladding  of  optical  fibers.  We  attempted,  without  success,  to  perform  absolute 
measurements.  When  we  calibrated  the  microscope  with  a  fiber  of  known  diameter, 
however,  we  could  measure  the  cladding  diameter  within  0.1  \xm. 

Key  words:  coherence;  gray  scale;  imaging;  micrometer;  microscope;  optical  fiber 
geometry;  optical  fibers;  video  microscopy 

1.    Introduction 

This  Technical  Note  describes  a  measurement  of  the  mean  diameter  and  noncircularity  of  the 
cladding  of  optical  fibers,  more  specifically,  an  attempt  to  perform  the  measurement  absolutely  with  a 
video  microscope.  Engineers  must  have  accurate  measurement  of  both  mean  diameter  and 
noncircularity  to  make  precise,  low-loss  connectors.  The  industry  requires  a  device  which  can  measure 
the  mean  cladding  diameter  with  an  accuracy  of  0.1  fxm.  The  video  microscope  is  the  most  common 
instrument  for  making  measurements  of  fiber  geometry. 

The  most  accurate  measurements  are  made  with  interferometric  microscopes  and  contact 
micrometers.  In  our  division,  we  use  a  contact  micrometer  to  make  diameter  measurements  accurate 
within  40  nm  [1].  We  manufacture  Standard  Reference  Material  (SRM)  fibers  for  the 
telecommunications  industry  so  that  they  have  calibration  artifacts  for  their  measuring  instruments.  The 
SRMs  are  approximately  125  jim  in  diameter  and  are  very  nearly  circular. 

We  measured  a  number  of  fibers  with  the  micrometer  and  used  them  to  check  our 
measurements,  we  did  not  succeed  in  making  the  video  microscope  an  absolute  device:  Our  mean 
diameter  measurements  were  consistently  high  by  0.19  ^m.  Nevertheless,  the  3a  uncertainty  in  our 
measurements  is  0.07  i^m,  which  is  better  than  the  required  accuracy.  Thus,  the  microscope  provides 
the  necessary  accuracy,  but  only  if  calibrated  with  an  artifact  such  as  the  SRM.  We  also  measured  the 
noncircularity  and  the  position  of  the  major  axis  of  an  ellipse  fitted  to  the  edge  of  the  fiber. 

The  video  microscope  has  a  number  of  inherent  problems.  As  with  all  optical  equipment,  the 
focus  of  the  instrument  is  critical.  Slight  defocusing  can  greatly  change  the  apparent  diameter  of  the 
optical  fiber  being  viewed.  Because  of  the  depth  of  focus  of  the  instrument,  it  is  difficult  to  find  the 
precise  focus.    Sensitivity  to  focus  also  creates  a  problem  when  the  fiber  end  face  is  not  perpendicular 
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to  the  optical  axis  of  the  microscope.  When  the  end  is  tihed,  parts  of  the  fiber  are  in  focus  while  other 
parts  are  not  in  focus.  This  creates  the  illusion  of  a  noncircular  fiber  end.  Even  when  the  image  is 
correctly  focussed,  the  image  is  not  sharp  because  of  diffraction  or  aberrations. 

The  width  of  the  image  of  the  edge  is  on  the  order  of  1  fim  as  a  result  of  diffraction  by  the 
aperture  stop  of  the  objective.  To  find  the  intensity  which  corresponds  to  the  edge  of  the  fiber,  we  must 
use  the  theory  of  partially  coherent  light.  The  width  of  the  edge  profile  and  the  circular  nature  of  the 
fiber  end  combine  with  the  rectangular  geometry  of  the  pixels  in  the  ccd  camera  to  create  a  systematic 
deviation  of  the  apparent  edge.  Precise  characterization  of  the  deviation  could  result  in  a  more  accurate 
location  of  the  edge  of  the  fiber. 


2.   Apparatus 

We  inherited  the  microscope  used  by  Steven  Mechels  [2]  and  Lakshmikant  Mamileti  [3]. 
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Figure  1. 


Detail  of  illuminator  and  microscope  objective. 


The  microscope  itself  consists  of  an  objective  lens  and  a  ccd  camera  connected  by  a  metal  tube. 
Inside  the  metal  tube  are  two  glare  stops.  Four  hardened  steel  rods  rigidly  hold  the  objective,  tube,  and 
camera.    The  objective  has  a  numerical  aperture  of  0.65  and  a  magnification  of  40x.   We  chose  the  40x 


magnification  to  ensure  that  the  image  of  the  fiber  would  fill  about  80%  of  the  available  pixel  array  in 
the  camera. 

The  objective  is  designed  to  function  with  a  coverslip.  Use  of  the  objective  without  a  coverslip 
results  in  spherical  aberration  [4],  so  we  attached  a  coverslip  to  a  cylinder  of  phenolic.  The  cylinder 
fits  over  the  object  and  inside  the  ring  illuminator  as  shown  in  Figure  1. 

The  camera  is  a  frame  transfer  array  as  discussed  in  [5].  It  has  512  pixels  horizontally  and  220 
vertically.  The  pixels  are  rectangular  and  measure  1 1  jim  horizontally  and  27  fxm  vertically.  Because 
of  the  way  the  microscope  is  calibrated,  the  absolute  dimensions  of  the  pixels  are  not  important. 

The  camera  reads  the  information  in  the  pixels  one  line  at  a  time.  A  set  of  220  lines  is  called 
a  field.  Our  camera  operates  in  the  interlaced  mode:  it  reads  two  separate  fields.  One  field  consists 
of  lines  corresponding  to  the  pixels  themselves.  The  second  field  consists  of  lines  of  pseudopixels. 
Each  pseudopixel  is  centered  vertically  between  the  real  pixel  above  and  the  real  pixel  below.  The 
light-gathering  area  of  the  pseudopixel  is  half  of  the  pixel  above  and  half  of  the  pixel  below.  The 
camera  achieves  this  second  field  by  biasing  the  pixels.  The  net  image  is  a  512  x  440  array  called  a 
frame. 

The  output  of  the  camera  goes  to  a  commercially  available  frame  digitizer  running  on  a  personal 
computer  and  to  a  video  monitor.  The  resulting  array  assigns  values  between  0  and  255  for  each  pixel. 
We  use  the  monitor  to  visually  position  and  focus  the  image  of  the  fiber  end.  Light  from  a  ring 
illuminator  reflects  from  a  plate  behind  the  fiber  to  back-illuminate  it.  The  ring  illuminator  has  a  hole 
drilled  through  its  center  and  fits  snugly  around  the  objective.  Instead  of  locating  a  green  filter  between 
the  objective  and  the  camera  (that  is,  inside  the  microscope)  [2],  we  drilled  a  hole  in  a  green  filter 
which  passes  a  band  of  light  centered  at  525  nm  and  attached  it  to  the  ring  illuminator.  We  use  a 
narrow  band  to  reduce  chromatic  aberration.  Attaching  the  filter  to  the  ring  illuminator  allows  us  to 
illuminate  the  core  with  light  of  any  wavelength. 

The  fiber  is  held  in  place  by  a  three-axis  optical  stage.  The  stage  has  micrometers  to  control 
position  along  the  optical  axis  of  the  microscope  as  well  as  two  orthogonal  directions  perpendicular  to 
the  optical  axis.  We  will  call  the  axis  parallel  to  the  optical  axis  the  z  axis  and  the  perpendicular 
directions  the  x  and  y  axes.  The  x  axis  is  parallel  to  the  floor.  The  perpendicular  controls  position  the 
fiber  in  the  center  of  the  field  of  view  of  the  microscope.  The  parallel  controls  focus  the  image  of  the 
fiber.  There  are  two  micrometers  for  focus  control.  One  micrometer  is  marked  in  25  \xm  (0.001  in) 
increments,  the  second  in  1  \xm  (40  fiin)  increments. 

To  accommodate  an  SRM  fiber,  we  attached  a  machined  aluminum  box  to  the  back  of  the  white 
aluminum  plate.  The  SRM  is  secured  in  the  box  by  two  screws.  The  fiber  can  be  extended  through 
a  hole  in  the  aluminum  plate.  The  plate  reflects  light  as  a  source  of  back-illumination.  The  plate  is 
fixed  to  a  beam  splitter  mount  which  is  attached  to  the  optical  stage,  as  shown  in  Figure  2. 
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Figure  2. 


Beam  splitter  mount  and  white  plate. 


The  beam  splitter  mount  provides  control  over  the  relative  angle  between  the  plane  of  the  fiber 
end  and  the  optical  axis.  The  two  possible  rotations  are  about  the  x  axis  (pitch)  and  the  y  axis  (yaw). 
We  modified  the  system  by  attaching  a  black  mask  to  the  aluminum  plate.  To  fix  the  size  of  the 
illumination  spot,  we  punched  a  0.95  cm  diameter  hole  in  the  center  of  the  mask.  The  hole  in  the  mask 
is  concentric  with  the  hole  in  the  metal  plate.  We  punched  the  hole  in  the  mask  both  to  fix  the  size  of 
the  illumination  spot  and  to  make  it  independent  of  the  position  of  the  ring  illuminator. 

An  SRM  fiber  is  enclosed  in  an  aluminum  box,  or  SRM  holder,  to  protect  it.  The  fiber  is  glued 
into  a  ferrule  in  a  brass  barrel;  the  barrel  itself  rests  in  the  aluminum  box.  The  front  of  the  box  has  a 
hole  though  which  the  fiber  projects.  The  barrel  slides  forward  and  backward  to  project  and  retract  the 
fiber  end.  The  barrel  is  marked  in  45°  increments  and  rotates  inside  the  box.  All  of  the  fibers  were 
mounted  in  SRM  holders. 


We  also  illuminate  the  core  of  the  fiber  to  investigate  its  geometry.  We  use  a  tungsten  lamp, 
a  lens,  and  a  second  microscope  objective  to  focus  light  into  the  core.  We  wrap  the  fiber  around  a 
mandrel  to  scramble  the  modes.  The  fibers  used  are  single-mode  fibers  for  1.3  }xm  light;  they  are 
muitimode  at  visible  wavelengths.  Multimode  operation  gives  rise  to  a  nonuniform  intensity  profile  in 
the  core;  this  can  create  an  error  in  measuring  the  positions  of  the  edge  and  the  center  of  the  core. 
Wrapping  the  fiber  about  a  mandrel  scrambles  the  modes  and  causes  a  more  uniform  intensity  profile. 


Figure  9  is  the  angle  of  the  major  axis  as  a  function  of  tilt.    Yaw  draws  the  major  axis  toward 
■90°  (vertical)  whereas  pitch  draws  it  toward  0°  and  —180°  (horizontal). 
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Figure  9 

M 

ajor 

axis  angle 

as  a  function  of  tilt  angle. 

For  the  second  set  of  data,  we  used  a  more  severely  noncircular  fiber  (Figures  10,  11,  and  12). 
The  contact  micrometer  measured  the  mean  diameter  as  125.522  |am  ±  0.04  fim,  the  noncircularit>  as 
0.25%,  and  the  angle  as  —89°.    These  data  also  support  the  suppositions. 
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The  graph  of  mean  diameter  as  a  function  of  tilt  angle  in  Figure  10  is  less  smooth  than  that  in 
Figure  7.  The  point  labeled  0°  yaw  on  the  graph  is  no  longer  0°  yaw,  possibly  because  the  yaw  screw 
was  very  loose  and  near  the  end  of  its  travel.  As  in  Figure  7,  an  increment  of  yaw  has  greater  effect 
than  an  equal  increment  of  pitch. 
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Figure  10.         Mean  diameter  as  a  function  of  tilt  angle. 
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The  noncircularity  in  Figure  1 1  has  a  minimum  at  0°  yaw.  It  has  a  local  maximum  at  0°  pitch 
and  local  minima  at  9°  and  -6°.  Just  as  for  Figure  8,  these  are  the  results  we  expect,  since  the  major 
axis  of  the  fiber  is  near  -90°.  At  the  local  minima,  the  vertical  diameter  and  the  horizontal  diameter 
are  the  same;  this  places  the  major  axis  near  —45°  or  -135°. 
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Figure  12  shows  that  the  angle  of  the  major  axis  is  —60°  at  9°  pitch  and  -120°  at  -6°  pitch. 
These  are  close  to  the  predictions  of  —45°  and  -  135°.  For  larger  pitch  angle,  the  major  axis  is  drawn 
to  the  horizontal.   Yaw  does  not  affect  the  direction  of  the  major  axis  since  it  was  already  nearly  —90°. 
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Figure  12.         Angle  of  major  axis  as  a  function  of  tilt  angle. 


The  effects  of  tilting  the  fiber  do  not  change  as  the  fiber  rotates  about  its  axis.  This  allows  us 
minimize  the  tilt.  By  looking  at  an  edge  table  in  polar  coordinates,  we  can  determine  whether  the  angle 
of  the  major  axis  rotates  as  the  fiber  rotates. 

In  a  plot  of  radius  as  a  function  of  angle,  a  circular  fiber  appears  as  a  horizontal  line  because 
the  radius  is  independent  of  angle.  An  elliptical  fiber  appears  as  a  sinusoid.  The  two  maxima 
correspond  to  the  major  axis;  the  two  minima  correspond  to  the  minor  axis.  Figure  13  is  a  plot  of 
radius  as  a  function  of  angle  for  a  very  nearly  circular  fiber. 
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Figure  13.         Radius  of  fiber  as  a  function  of  angle. 

Ideally,  if  the  fiber  rotates  by  45°,  then  the  maxima  and  minima  should  also  shift  by  45°.  The 
magnitudes  of  the  maxima  and  minima  are  determined  by  the  noncircularity  of  the  fiber.  Tilting  causes 
a  second  sinusoid  to  be  added  to  the  noncircularity  sinusoid.  The  maxima  and  minima  from  the  tilt 
sinusoid  do  not  rotate  with  the  fiber. 

If  there  is  pitch,  the  maxima  for  the  second  sinusoid  appear  at  0°  and  1 80°.  If  there  is  yaw,  the 
maxima  are  at  90°  and  270°.  If,  as  the  fiber  rotates,  there  are  persistent  maxima  at  90°  and  270°.  yaw 
must  be  adjusted.  If  the  persistent  maxima  appear  at  0°  and  180°,  pitch  must  be  adjusted.  By  checking 
the  movement  of  the  maxima,  we  can  adjust  pitch  and  yaw  to  0°  within  1°  or  so.  Inside  this  range,  we 
can  no  longer  distinguish  between  two  edge  tables. 

It  is  difficult  to  determine  the  precise  effect  of  tilt  on  the  mean  diameter  (Figures  7  and  10). 
The  uncertainty  due  to  defocus  is  so  large  that  it  overwhelms  the  effect  of  tilting.  Therefore,  we  fitted 
a  parabola  to  the  mean  diameter  as  a  ftinction  of  tilt  angle.  The  value  of  the  fitted  curve  drops  by 
0.004  |j,m  when  the  tilt  angle  is  1°.  We  will  assume  that  the  best  value  is  at  the  middle  of  this  range 
and  that  the  distribution  about  the  best  value  is  Gaussian.  If  we  assume  the  0.004  f^m  range  to  be  6a 
then  the  la  uncertainty  is  ±0.0007  |im.  The  range  of  the  noncircularity  is  about  0.04%  over  a  2°  range; 
again  we  will  assume  that  this  range  is  6a.  The  la  uncertainty  of  the  noncircularity  due  to  tilt  is  thus 
±0.007%.    The  range  in  the  angle  of  the  major  axis  is  8°;  the  la  uncertainty  is  ±1.3°. 
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However,  the  beam  splitter  mount  is  not  the  only  possible  source  of  tilt.  Cleaving  of  the  fiber 
end  does  not  always  leave  a  surface  perpendicular  to  the  axis  of  the  fiber.  When  we  cleave  the  fibers, 
we  check  the  flatness  and  perpendicularity  of  the  ends  with  a  commercial  Michelson  interferometer  that 
uses  730  nm  light.  The  best  ends  display  about  three  fringes;  this  implies  an  angle  of  about  0.5°,  which 
is  enough  to  alter  the  appearance  of  the  fiber  in  the  video  microscope.  It  is  impossible  to  tell  the 
relative  geometry  between  the  tilt  of  the  cleave  and  the  actual  major  axis  of  the  fiber.  The  tilt  axis 
could  be  parallel  to  the  major  axis  of  the  fiber,  perpendicular  to  the  major  axis,  or  somewhere  in 
between.  Because  this  tilt  is  part  of  the  fiber  end,  we  cannot  adjust  for  it.  The  best  we  can  do  is  to 
try  to  align  the  fiber  with  the  optical  axis.  The  uncertainty  due  to  the  cleave  angle  is  about  the  same 
as  that  due  to  the  tilt  of  the  beam  splitter  mount,  since  we  cannot  determine  the  geometry  of  the  tilt 
angle  and  the  true  major  axis. 


6.   Partial  Coherence 

Even  a  perfectly  focused  image  is  not  perfectly  sharp  due  to  diffraction  by  the  aperture  stop  of 
the  objective;  this  causes  the  image  of  the  edge  to  have  a  finite  width.  The  camera  records  the  intensity 
pattern  of  the  edge  in  the  image  plane  of  the  microscope  which  we  call  the  edge  profile.  The  width  of 
the  edge  profile  makes  it  unclear  where  the  true  edge  lies.  To  calculate  the  relative  intensity  at  the 
edge,  we  used  the  theory  of  partially  coherent  light  discussed  by  Beran  and  Parrent  [8],  Bom  and  Wolf 
[9],  and  Goodman  [10]. 

The  ccd  detects  the  time-averaged  intensity;  intensity  is  the  square  of  the  complex 
electromagnetic  field.  Consider  two  point  sources  of  radiation.  The  total  field  at  any  point  in  space 
is  the  sum  of  the  fields  due  to  the  two  points.   The  ccd  sees  the  time  average  of  the  square  of  the  sum, 

L  -  <\E,  +  EJ^>  =  <\E,\^>  +  <\EJ^>  +  <£,£  +  E,E,>.  (1) 

time  avg  II  21  111  i2l  12  12 

When  the  square  is  expanded,  the  first  two  terms  are  then  intensities  due  to  the  individual  sources. 
However,  there  is  also  a  cross  term.  If  there  is  a  fixed  phase  relationship  between  the  fields  due  to  the 
two  points,  then  the  time-averaged  cross  term  can  have  a  value  other  than  0.  The  cross  term  will  be 
a  maximum  if  the  fields  are  in  phase  and  0  if  they  are  out  of  phase.  The  light  is  called  coherent,  and 
we  must  work  with  the  electromagnetic  field.  If,  on  the  other  hand,  there  is  no  correlation  between  the 
phases  of  the  two  fields,  then  the  time  average  causes  the  cross  term  to  vanish.  This  is  called  the 
incoherent  case,  and  we  add  the  individual  intensities  to  get  the  final  intensity. 

Unfortunately,  for  extended  sources  imaged  through  finite  apertures,  the  light  is  neither  coherent 
nor  incoherent  but  is  partially  coherent.  The  convenient  quantity  is  neither  the  electromagnetic  field 
nor  the  intensity  but  the  mutual  coherence  fiinction 


T{x,y;w,z;r)  -  <Vix,y,t+x)Viw,z,t)>. 


(2) 


Partially  coherent  light  is  linear  in  the  mutual  coherence  function.  V  is  the  complex  scalar  field  at  a 
point  (x,y).  The  point  (w,z)  is  in  the  same  plane  as  the  point  (x,y).  The  bar  indicates  the  complex 
conjugate,  x  is  a  time  difference.  The  o  indicate  a  time  average  in  t.  The  intensity  is  the  mutual 
coherence  function  at  a  single  point  with  x=0;  that  is, 

I{x,y)  -  nx,y;x,y;0).  (3) 
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The  mutual  coherence  function  satisfies  two  wave  equations,  one  for  each  set  of  spatial 
coordinates.  Therefore,  if  we  know  the  mutual  coherence  function  in  the  object  plane  and  the  Green's 
function  (transfer  function)  of  the  optical  system,  we  can  find  the  mutual  coherence  function  in  the 
image  plane.  The  Green's  function  for  the  mutual  coherence  function  in  one  set  of  spatial  variables  is 
the  same  as  the  Green's  function  for  the  electromagnetic  field  under  coherent  illumination.  With  the 
mutual  coherence  function  in  the  image  plane,  we  can  find  the  intensity.  Both  coherent  and  incoherent 
illumination  can  be  described  as  certain  limits  of  partially  coherent  illumination. 

The  mutual  coherence  function  of  a  circular,  mutually  incoherent  source  can  be  described  by 
a  characteristic  distance.  A  mutually  incoherent  source  is  one  in  which  there  is  no  correlation  between 
the  phases  of  the  fields  arising  from  any  two  points.  In  other  words,  the  mutual  coherence  function  of 
the  source  is  a  delta  fiinction  of  position.  The  characteristic  distance  is  the  distance  between  the 
maximum  and  first  zero  of  the  mutual  coherence  function.  This  distance,  the  radius  of  coherence,  is 
determined  by  the  numerical  aperture  of  the  source.  In  the  paraxial  approximation,  the  radius  of 
coherence  is 

r    =  0.61—,  (4) 

NA  ^  ' 

s 

where  NA^  is  the  numerical  aperture  of  the  source.  NA^  is  the  sine  of  the  angle  formed  by  the  optical 
axis  and  a  line  from  the  fiber  end  to  the  edge  of  the  source.  In  Figure  14,  NA^  =  sin(a).  The  phases 
of  points  inside  this  distance  are  highly  correlated. 
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Figure  14.         Numerical  apertures. 

The  resolution  limit  of  the  microscope  is  a  second  characteristic  distance  which  influences  the 
degree  of  coherence.  The  resolution  limit  is  the  minimum  distance  between  two  equally  intense  point 
sources  that  appear  as  separate  points  in  the  image.  The  resolution  limit  in  the  object  plane  is  defined 
as 


BL  =  0.61 


NA, 


(5) 


where  NA,,  is  the  numerical  aperture  of  the  objective.    In  Figure  14,  NA^  =  sin(b). 

The  degree  of  coherence  of  the  system  is  determined  by  these  two  distances.  If  the  coherence 
length  is  much  larger  than  the  resolution  limit,  then  the  system  appears  to  be  highly  coherent.  If  the 
coherence  length  is  much  smaller  than  the  resolution  limit,  then  the  system  appears  incoherent. 

We  define  the  coherence  parameter  a  [11]  as  the  ratio  of  the  resolution  limit  to  the  coherence 
length  or,  equivalently,  the  ratio  of  the  numerical  aperture  of  the  source  to  the  numerical  aperture  of 
the  microscope  objective.  A  coherent  source  is  a  point  source  and  has  a  coherence  parameter  of  0.  By 
contrast,  a  source  can  be  considered  incoherent  only  if  the  objective  has  a  numerical  aperture  of  0.   An 
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We  use  white  light  to  illuminate  the  core  to  increase  the  number  of  modes;  this  also  causes  a  more 
uniform  intensity  profile  in  the  core.  We  cannot  illuminate  the  core  with  a  single  symmetric  mode 
because  the  ccd  is  not  sensitive  to  1 .3  |j.m  light. 

All  of  the  apparatus  except  the  computer  and  the  video  monitor  is  held  in  line  by  an  optical 
bench. 


3.   Calibration  and  Data  Reduction 

To  measure  diameter  accurately,  the  optical  system  must  be  calibrated.  That  is,  we  need  to 
determine  the  correspondence  between  numbers  of  pixels  in  the  camera  and  distance  in  the  object  plane. 
We  calibrated  the  microscope  by  using  a  chromium-on-glass  reticle.  The  reticle  has  21  clear  lines  on 
an  opaque  background.  The  lines  are  approximately  10  \xm  apart.  The  reticle  had  been  sent  previously 
to  NIST  Gaithersburg  for  accurate  measurement  of  the  line  separations. 

To  calibrate  horizontally,  we  digitize  an  image  with  the  lines  running  vertically  on  the  screen. 
They  are  separated  by  spaces  in  the  x  direction.  The  image  has  15  peaks  in  intensity  [2].  We  calculate 
the  distance  in  pixels  between  the  lines  by  calculating  the  distance  between  the  peaks.  Because  each 
of  the  peaks  has  a  finite  width  and  is  nearly  symmetric  about  its  center,  we  find  the  center  of  the  line 
by  finding  the  center  of  the  segment  between  the  half-maximum  intensity  points.  The  distance  between 
two  lines  is  the  distance  between  their  centers.  By  fitting  a  polynomial  in  pixel  number  to  the  data,  we 
get  distance  as  a  function  of  pixels.  Only  the  linear  term  is  significant;  it  is  a  scale  factor  in 
micrometers  per  pixel.  We  repeat  the  same  procedure  with  the  lines  rotated  by  90°  for  the  vertical 
calibration.  The  calibration  is  described  in  more  detail  in  [2].  The  second  term  of  the  polynomial  fit 
discussed  in  that  reference  turns  out  to  be  unreproduceable;  its  inclusion  results  in  a  distorted  fiber 
image. 

The  final  output  of  this  system,  a  512  x  440  array  of  hexadecimal  numbers,  represents  the  image 
of  the  fiber  end  in  256  gray  levels.  Jack  Wang  and  Dominic  Vecchia  of  the  Statistical  Engineering 
Division  at  NIST  wrote  the  software  which  locates  the  edge  and  fits  an  ellipse  to  the  edge  [3]. 

The  output  of  the  program  contains  the  Cartesian  coordinates  of  the  points  on  the  edge,  the 
coordinates  of  the  center  of  the  ellipse,  the  semi-major  axis,  the  semi-minor  axis,  and  the  angle  between 
the  semi-major  axis  and  the  horizontal  axis,  as  well  as  the  same  quantities  for  the  core  of  the  fiber. 

The  Cartesian  coordinates  of  the  points  on  the  edge  are  called  an  edge  table.  There  is  an  edge 
table  for  the  cladding  and  one  for  the  core. 

To  get  the  x,y  data  from  the  array,  the  program  splits  the  data  into  two  sections.  The  first 
section  includes  the  fiber  edge  from  -45°  to  45°  and  135°  to  225°.  In  this  section,  the  data  are  scanned 
horizontally  line  by  line.  The  edge-finding  routine  is  flexible  and  can  find  the  edge  at  a  specific 
fraction  of  the  maximum  intensity  or  at  the  inflection  point.  The  second  section  includes  the  fiber  edge 
from  45°  to  135°  and  225°  to  315°,  and  the  program  scans  the  data  in  vertical  lines. 

The  program  fits  an  ellipse  to  the  points  in  the  edge  table  in  two  stages.  The  program  first  uses 
least  median  of  squares  (LMS)  to  remove  outliers  and  then  least  sum  of  squares  (LSS)  to  fit  an  ellipse 
to  the  remaining  points.    This  procedures  differs  from  the  procedure  used  b>  [3]. 


LMS  is  a  robust  method  that  minimizes  the  median  of  the  squares  of  the  residuals  and  is 
insensitive  to  outliers  [6].  A  residual  is  the  distance  between  a  data  point  and  the  fitted  ellipse,  though 
that  distance  need  not  correspond  to  a  physical  distance,  depending  on  the  fitting  algorithm. 

The  program  offers  three  options  for  the  LSS  fit;  each  uses  a  different  definition  of  the  distance. 
The  first  option  is  called  "convenient."  It  uses  a  residual  which  is  convenient  for  calculation.  This 
residual  is  not  the  Euclidean  distance  between  the  point  and  the  ellipse.  The  second  option,  called 
"perpendicular,"  uses  the  perpendicular  Euclidean  distance  between  the  point  and  the  ellipse.  The  third 
option  is  called  "rigorous."  Here  we  assume  that  the  horizontal  (x)  coordinate  of  an  edge  point  from 
a  vertically  scanned  line  of  data  has  no  uncertainty  and  define  the  residuals  in  the  vertical  (y)  direction. 
Likewise,  we  define  the  residuals  of  a  horizontally  scanned  line  in  the  vertical  direction.  We  had  hoped 
that  the  "rigorous"  fits  would  substantially  improve  the  fit.  However,  for  the  images  we  used,  the  noise 
was  small,  and  the  fits  were  not  significantly  different.  The  difference  between  the  convenient  LSS  and 
the  perpendicular  Euclidean  distance  LSS  is  discussed  more  rigorously  in  [7]. 

The  LMS  stage  does  not  radically  improve  the  fit  to  an  undamaged  fiber  end.  However,  if  the 
fiber  end  is  damaged,  the  LMS  stage  removes  the  points  at  the  damaged  part  of  the  fiber  and  still 
achieves  a  good  fit  [3].  The  ends  of  the  fibers  we  used  had  no  noticeable  defects  in  the  edge;  the  fits 
from  the  three  methods  were  not  significantly  different. 

In  addition  to  the  edge  table  in  Cartesian  coordinates,  we  found  it  useful  to  look  at  the  edge  in 
polar  coordinates.  In  a  plot  of  fiber  radius  as  a  function  of  angle  (an  r-0  plot),  an  ellipse  appears  as  a 
sinusoid  with  two  maxima  and  two  minima.  The  maxima  are  the  radii  along  the  major  axis  of  the 
ellipse.  Noncircularity  is  easier  to  detect  in  a  plot  of  this  type.  The  scale  can  be  expanded  about  the 
average  radius  to  show  fine  detail  in  the  edge  table. 


4.   Defocus 

Our  ability  to  focus  the  microscope  was  the  largest  source  of  uncertainty  in  the  diameter 
measurements.  The  image  on  the  monitor  appeared  sharp  over  a  focusing  range  of  2  fim.  We  call  this 
range  the  operating  depth  of  focus  to  distinguish  it  from  the  theoretical  depth  of  focus  of  the  objective 
alone.  We  looked  at  the  mean  diameter  of  a  fiber  end  as  a  function  of  defocus.  Define  positive 
defocus  when  the  fiber  is  too  close  to  the  objective,  negative  defocus  when  the  fiber  is  too  far  from  the 
objective.  Best  focus  is  the  middle  of  the  operating  depth  of  focus.  To  find  the  best  focus,  we  start 
with  the  fiber  out  of  focus  and  move  it  toward  focus  until  the  image  stops  changing.  We  call  this  point 
the  edge  of  focus.  If  we  approach  from  too  far  away,  we  find  the  back  edge  of  focus.  If  we  approach 
from  too  near,  we  find  the  forward  edge  of  focus.  We  then  move  the  fiber  1  fj.m,  one  half  the  operating 
depth  of  focus,  into  focus.  The  fiber  is  now  near  the  middle  of  the  operating  depth  of  focus.  Figure 
3  is  a  graph  of  mean  diameter  as  a  function  of  defocus  for  two  separate  runs. 
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Figure  3. 


Mean  diameter  as  a  function  of  defocus. 


In  the  first  set  of  data  for  mean  diameter  as  a  function  of  defocus,  we  moved  the  fiber  back  to 
the  best  focus  before  taking  it  out  of  focus  for  each  point  (Run  1 ).  In  other  words,  we  found  the  best 
focus  and  digitized  the  image.  We  moved  the  fiber  1  ^im  out  of  focus  and  digitized  another  image. 
Then,  we  moved  the  fiber  to  best  focus,  moved  it  2  |u.m  out  of  focus,  and  digitized  a  third  image,  and 
so  on.    To  eliminate  backlash,  we  always  found  the  best  focus  by  first  finding  the  back  edge  of  focus. 

In  Run  2,  we  first  found  the  back  edge  of  focus  and  then  best  focus.  We  moved  the  fiber 
forward,  without  returning  to  focus  as  in  Run  1 ,  to  acquire  positive  defocus  data.  We  then  found  the 
forward  edge  of  focus  and  best  focus.  We  moved  the  fiber  away  from  the  objective,  without  returning 
to  focus,  to  acquire  negative  defocus  data.  Images  that  had  net  defocus  greater  than  2  |im  were  visibK 
out  of  focus. 


Both  sets  of  data  have  very  similar  shapes.  The  best  focus  appeared  to  be  at  a  local  maximum 
of  the  mean  diameter.  However,  we  find  the  best  focus  in  a  subjective  manner  and  cannot  resolve  the 
curve  accurately  enough  to  state  positively  that  the  best  focus  is  a  local  maximum.  The  mean  diameter 
has  a  range  of  0.3  )im  over  the  operating  depth  of  focus.  However,  0.3  fim  is  an  o\erestimation  of  the 
uncertainty  due  to  focussing.    We  can  repeatably  find  the  best  focus  well  within  the  operating  depth  of 


focus.  We  made  12  runs  in  order  to  find  a  statistical  estimate  of  the  uncertainty.  For  each  run,  we 
found  the  best  focus  starting  at  the  back  edge  and  starting  at  the  front  edge.  Figure  4  represents  the 
mean  diameter  for  each  of  the  12  runs.  The  mean  diameter  of  all  of  the  measurement  is  125.122  fim. 
The  overall  spread  of  values  is  0.1  fim  counting  an  apparent  outlier  at  Run  7.  Figure  3  shows  a 
maximum  in  the  mean  diameter  at  0  defocus,  so  we  will  assume  that  the  region  near  best  focus  is  nearly 
flat.  We  will  also  assume  that  the  distribution  of  values  about  the  mean  is  Gaussian.  Without  Run  7, 
the  standard  deviation  or  la  uncertainty  of  the  mean  diameter  is  ±0.018  |im.  The  3a  uncertainty  is 
±0.054  fim,  well  within  the  goal  of  0.1  fxm. 
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Figure  4. 


Mean  diameter  at  best  focus. 


The  mean  diameter  of  the  fiber  is  124.935  \x.m  with  a  3a  uncertainty  of  ±0.04  |j,m  as  measured 
on  the  contact  micrometer.  The  video  microscope  therefore  displays  an  error  of  0.19  fa.m,  which 
exceeds  the  goal  of  0.1  jim  absolute  accuracy.  The  la  uncertainty  in  this  error  is  the  standard  deviation 
of  the  mean  of  the  data,  or  ±0.004  |im.    The  3a  uncertainty  of  the  error  is  ±0.012  \x,m. 


Figure  5  shows  the  noncircularity  for  each  of  the  runs.  The  average  noncircularity  was  0.054%. 
Neglecting  Run  7  and  assuming  that  the  distribution  about  the  mean  is  Gaussian,  the  la  uncertainty  is 
±0.003%. 
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Figure  5.  Noncircularity  at  best  focus. 


Figure  6  shows  the  angle  between  the  major  axis  and  the  horizontal  plane  for  each  run.  The 
average  angle  is  —80°  and,  by  reasoning  similar  to  that  for  the  mean  diameter,  the  Ict  uncertainty  is 
±8°. 
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Figure  6.  Angle  of  major  axis  at  best  focus. 


5.   Tilt 

When  the  fiber  end  is  not  perpendicular  to  the  optical  axis,  the  data  of  Figure  3  suggest  that 
the  out-of-focus  sections  of  the  end  will  appear  too  small,  whereas  the  in-focus  sections  should  remain 
unchanged.  The  net  effect  of  the  tilt,  therefore,  is  a  smaller  mean  diameter,  an  altered  noncircularity, 
and  a  shift  in  the  apparent  position  of  the  major  axis. 

Consider  a  fiber  end  tilted  about  some  axis.  The  diameter  parallel  to  that  axis  will  appear 
unchanged,  while  the  diameter  perpendicular  to  the  axis  will  appear  too  small.  As  the  tilt  angle 
becomes  larger,  the  apparent  difference  between  the  two  diameters  will  increase.  Regardless  of  the  true 
position  of  the  major  axis,  the  apparent  major  axis  will  ultimately  align  with  the  tilt  axis.  The 
noncircularity  can  follow  either  of  two  paths.  If  the  diameter  parallel  to  the  tilt  axis  is  the  major  axis 
of  the  fiber  end  face,  the  noncircularity  will  increase  with  increasing  tilt.  If,  however,  that  diameter  is 
the  minor  axis,  the  noncircularity  will  first  decrease  as  the  apparent  perpendicular  diameter  is  brought 
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nearer  the  value  of  the  parallel  diameter.    Then,  once  the  perpendicular  diameter  is  made  smaller  than 
the  parallel  diameter,  the  noncircularity  will  increase  once  again. 

We  tested  these  suppositions  by  using  the  angle  adjustments  on  the  beam  splitter  mount.  The 
beam  splitter  mount  has  two  screws  which  control  two  orthogonal  rotations.  The  distance  between  the 
screws  and  their  pivots  is  32.8  mm  (1.29  in).  The  screws  have  80  threads  per  inch.  Let  us  say  that 
positive  pitch  brings  the  top  edge  of  the  fiber  away  from  the  objective  and  that  positive  yaw  brings  the 
right  edge  of  the  fiber  end  away  from  the  objective.  When  we  took  data  changing  yaw,  we  focused  on 
the  top  and  bottom  of  the  fiber  image  because  these  do  not  change  their  relative  positions  when  we 
change  the  yaw  angle.  When  we  changed  pitch,  we  focused  on  the  left  and  right  edges.  We  always 
came  from  the  back  edge  of  focus. 

In  the  first  set  of  data,  we  used  a  very  nearly  circular  fiber,  whose  mean  diameter  is 
124.935  fim,  noncircularity  is  0.3  %,  and  major  axis  angle  is  —70°.  Our  predictions  are  all  supported 
by  the  data. 

Figure  7  is  the  mean  diameter  as  a  function  of  tilt.  As  with  defocus  in  Figure  3,  the  mean 
diameter  appears  to  have  a  maximum  at  0°  tilt.  For  the  same  angle,  yaw  appears  to  have  a  greater 
effect  on  the  mean  diameter  than  pitch.  This  could  be  due  to  the  different  pixel  geometries  in  the  two 
directions.    (The  travel  of  the  screw  limited  the  amount  of  negative  yaw.) 
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Figure  7.  Mean  diameter  as  a  ftinction  of  tilt  angle. 
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Figure  8  is  the  noncircularity  as  a  function  of  tilt  angle.  The  noncircularity  appears  to  have  a 
minimum  at  0°  yaw.  Because  the  major  axis  of  the  fiber  is  near  —90°,  we  expect  the  noncircularity 
to  decrease  at  small  pitch;  in  fact,  it  decreases  slightly  at  ±3°  pitch  and  then  increases.  Yaw  also 
appears  to  have  the  greater  effect  on  the  noncircularity. 
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Figure  8.  Noncircularity  as  a  function  of  tilt  angle. 
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to  a  period  of  two  pixels.   We  used  a  discrete  one-dimensional  convolution  to  simulate  its  effect  on  the 
image.    We  saw  no  effect  on  any  of  the  geometric  parameters  or  on  the  edge  table. 

We  also  did  a  simulation  to  see  whether  pixelating  changed  the  location  of  the  edge.  The  pixels 
are  on  average  spatial  low-pass  filters.  We  used  a  one-dimensional  square  convolution  kernel  to 
approximate  a  low-pass  filter  and  convolved  it  with  the  theoretical  edge  profile.  The  original  edge 
profile  and  the  shifted  edge  profiles  are  shown  in  Figure  22.  The  position  axis  is  expanded  about  0  to 
make  the  shift  in  the  intensity  at  the  edge  more  clear. 
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Figure  22.         Intensity  at  the  edge:  theoretical  and  pixelated 
1  resolution  limit  =  3.14. 

This  simulation  indicates  that  low-pass  filtering  by  the  pixels  increases  the  apparent  intensity 
at  the  edge  by  0.005.  The  effect  of  looking  at  an  intensity  higher  by  0.005  would  be  a  mean  diameter 
error  +0.01  \xm,  a  value  not  large  enough  to  account  for  the  apparent  discrepancy.  However,  since  the 
real  and  theoretical  edge  profiles  are  different,  this  simulation  is  meaningful  only  qualitatively. 


8.    Fiber  Rotation 

The  barrel  in  the  SRM  holder  is  free  to  rotate  and  is  marked  at  every  45°.    Rotating  the  fiber 
about  its  axis  is  usefial  to  differentiate  between  the  true  geometry  of  the  fiber  itself  and  artifacts  of  the 
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microscope.  Ideally,  the  noncircularity  and  the  mean  diameter  should  not  change  as  the  fiber  is  rotated. 
The  angle  of  the  major  axis  should  move  by  the  same  angle  as  the  fiber.  We  used  a  very  nearly 
circular  fiber,  noncircularity  0.03%,  to  investigate  the  effects  of  rotating  the  fiber. 

We  rotated  the  fiber  in  45°  increments  and  took  an  image  at  each  increment.  Figure  23  is  the 
angle  of  the  major  axis  as  a  function  of  the  angle  of  fiber  rotation.  The  angle  of  the  major  axis  moves 
with  fiber  rotation  but  not  as  we  anticipated.  When  the  rotation  angle  is  near  0°,  the  angle  of  the  major 
axis  is  very  close  to  what  was  measured  by  the  contact  micrometer.  However,  as  we  rotate  the  fiber, 
the  measured  angle  moves  away  from  our  expectation.  Initially,  we  thought  that  this  was  due  to  tilt  and 
adjusted  the  tilt  to  try  to  remove  the  glitch.  The  data  labeled  "center"  are  for  tilt  angle  as  close  to  0  as 
possible  using  the  calibration  method  described  before. 
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Figure  23.         Angle  of  major  axis  as  a  function  of  fiber  rotation  angle. 


As  Figure  23  shows,  we  could  minimize  the  glitch  but  not  remove  it.  We  think  it  comes  about 
because  the  axis  of  the  fiber  is  not  precisely  perpendicular  to  the  axis  of  rotation  of  the  barrel.  This 
is  verified  by  noting  that  the  image  of  the  end  moved  upward  slightly  at  0°  and  downward  at  180°. 

Because  the  fiber  is  not  parallel  to  the  optical  axis,  the  angle  between  the  fiber  end  and  the 
optical  axis  changes  as  the  fiber  is  rotated.  Since  we  minimize  tilt  with  the  fiber  at  0°,  the  tilt  will  not 
be  minimized  at  180°.  But  we  can  instead  minimize  the  average  tilt  for  all  rotations  by  minimizing  the 
glitch  in  the  curve  of  major  axis  angle  as  a  function  of  fiber  rotation  angle.    The  mean  diameter  as  a 


30 


function  of  the  fiber  rotation  angle,  Figure  24,  does  not  display  a  glitch  at  180°;  the  outlier  at  90°  is 
probably  a  blunder. 
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Figure  24.         Mean  diameter  as  a  function  of  fiber  rotation  angle. 

Figure  25  is  the  noncircularity  as  a  function  of  fiber  rotation  angle.  The  noncircularit\'  also 
changes  as  the  fiber  rotates.  There  appears  to  be  a  glitch  between  135°  and  225°.  However,  the 
uncertainties  due  to  the  cleave  and  tilt  are  large  enough  that  most  of  the  data  are  within  each  other's 
uncertainties. 
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Figure  25.         Noncircularity  as  a  function  of  fiber  rotation  angle. 


9.    Concentricity  Error 

Another  important  parameter  is  the  distance  between  the  center  of  the  cladding  and  the  center 
of  the  core.  This  is  also  known  as  the  concentricity  error.  We  looked  at  the  vector  between  the 
cladding  center  and  the  core  center  as  a  function  of  the  rotation  angle  of  the  fiber.  As  we  rotate  the 
fiber,  the  magnitude  of  the  vector  should  not  change,  but  its  angle  should  change  by  the  angle  of 
rotation.  That  is,  the  tip  of  each  vector  should  lie  on  a  circle  centered  about  the  center  of  the  cladding. 
Unfortunately,  the  magnitude  of  the  vector  changes  measurably,  and  the  angle  does  not  change  by  the 
expected  value.  Figures  26,  27,  and  28  show  the  tip  of  the  vector  that  starts  at  the  cladding  center  (0, 
0)  and  ends  at  the  core  center  for  three  sets  of  rotation  data. 
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Figure  26.  Cladding-core  vector  as  a  function  of  fiber  rotation  angle  (run  1). 
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Figure  27.  Cladding-core  vector  as  a  function  of  fiber  rotation  angle  (run  2). 
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Figure  28.  Cladding-core  vector  as  a  function  of  fiber  rotation  angle  (run  3). 

In  these  figures,  the  data  labels  represent  the  angle  of  rotation  of  the  fiber  with  respect  to  an 
arbitrary  starting  point.  The  origin  of  coordinates  is  the  center  of  the  cladding  as  calculated  from  an 
ellipse  fit.  In  each  case,  the  vector  rotates  only  very  approximately  in  a  circle,  but  the  center  of 
rotation  is  a  point  removed  from  the  origin  by  a  vector  whose  magnitude  is  approximately  0.26  fim  and 
whose  angle  is  —157°.    We  have  been  unable  to  identify  the  source  of  this  error. 


To  estimate  the  uncertainty  of  the  vector,  we  used  the  same  data  as  for  Figures  4,  5,  and  6. 
tips  of  the  concentricity  error  vectors  for  these  data  are  plotted  in  Figure  29. 


The 


34 


0.05  - 

o 

Oo 

o 
o 

o 

o 
o 
o 

oo         o 

o 
o 

-0.05- 

o 

o°    o 

o    o 
o 

)--0.15  - 

-0.25  - 

-0.35  - 

1 

1 

1 

1 

-0.40 

-0.30 

-0.20 

X 

-0.10 

0.00 

.Figure  29.    Core-cladding  vectors  for  data  sets  of  Figures  4,  5,  6. 


We  found  the  center  of  the  data  points  in  Figure  29  and  the  distance  between  the  center  and 
each  point.  The  standard  deviation  of  these  distances  is  an  estimator  of  the  uncertainty  of  the 
magnitude  of  the  vector;  it  is  equal  to  0.057  [xm.  The  largest  change  in  the  vector  angle  would  occur 
if  the  0.057  |im  deviation  were  perpendicular  to  the  offset  vector.  This  creates  a  la  angle  uncertainty 
of  12°. 

We  have  not  investigated  further,  but  presumably  to  calculate  the  concentricity  error  accurately, 
we  would  have  to  refer  the  measurement  to  the  tip  of  the  vector  rather  than  to  the  center  of  the 
cladding.    Something  along  these  lines  is  the  procedure  actually  used  in  commercial  instruments. 


10.    Summary  of  Uncertainties 

Table  1  contains  a  summary  of  the  uncertainties  and  their  sources.  These  uncertainties  are  for 
single  measurements  only.  The  artifact  uncertainty  comes  about  because  we  used  an  SRM  to  calibrate 
the  microscope  and  to  calculate  the  error  of  the  mean  diameter.  The  micrometer  measurement 
contributes  a  la  uncertainty  of  ±0.013  \xm  to  the  mean  diameter  measurement. 
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Table  1.    Summary  of  la  uncertainties. 


Mean 

diameter, 

|j,m 

Non- 
circularity, 
% 

Angle, 

0 

Vector 

magnitude, 

|j,m 

Vector 
angle,  ° 

Focus 

0.018 

0.003 

8 

0.057 

12 

Tilt 

0.0007 

0.007 

1.7 

Cleave  angle 

0.0007 

0.007 

1.7 

Edge-finding 
intensity 

0.0017 

0.000  07 

0.07 

Offset 

0.004 

Artifact 

0.013 

Total  (Ict) 

0.023 

0.01 

8.4 

0.057 

12 

Expanded 
(3a) 

0.07 

0.03 

25 

0.171 
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11.    Other  Fibers 

We  used  a  fiber  of  mean  diameter  124.935  fxm  and  noncircularity  0.03  %  to  calibrate  the 
microscope  and  investigate  the  uncertainties.  We  made  measurements  on  a  second  fiber  at  each  45° 
rotation  increment  to  check  the  results.  The  second  fiber  had  a  mean  diameter  of  125.522  |am,  a 
noncircularity  of  0.25  %,  and  a  major  axis  angle  of  —89°.  This  is  the  same  fiber  that  we  used  for  the 
second  run  of  tilt  data. 

Figure  30  is  a  graph  of  mean  diameter  as  a  function  of  fiber  rotation  angle.  The  mean  diameter 
has  been  adjusted  by  subtracting  the  0.19  |im  offset.  The  measured  values  are  consistently  below  the 
micrometer  measurement  but  mostly  within  the  3a  uncertainties  of  our  measurements. 
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incoherent  source  has  an  infinite  coherence  parameter.  The  mutual  coherence  function  of  a  nearly 
monochromatic,  circular  source  in  a  plane  a  finite  distance  away  can  be  found  using  the  van  Cittert- 
Zemike  theorem. 

We  used  a  half-plane,  opaque  for  y  <  0,  to  simulate  the  edge  of  the  fiber  (in  our  system,  the 
fiber  is  back-lighted  and  appears  as  a  dark  circle).  In  order  to  simplify  interpretation,  we  normalized 
all  spatial  variables  by  expressing  them  in  units  of  the  resolution  limit  of  the  microscope  in  the  object 
plane.  The  resolution  limit  in  the  image  plane  is  the  resolution  limit  in  the  object  plane  times  the 
magnification  of  the  microscope.  Two  objects  separated  by  one  resolution  limit  in  the  object  plane  will 
be  separated  by  one  resolution  limit  in  the  image  plane.  This  automatically  takes  into  account  the 
magnification  of  the  microscope.  It  also  makes  the  coherence  parameter  an  explicit  variable  in  the 
integral.  We  also  normalized  the  intensity  by  dividing  it  by  the  intensity  at  the  origin  in  the  absence 
of  the  half-plane.    The  intensity  in  the  image  plane  is 

I(a,b)  =  C  \dx\dy\dw\dz  somb{oir^somb{r^somb{r^,  (^) 


where  C  is  a  normalization  constant  and  the  sombrero  function  [12]  is 

2^,(^) 
somber)  =  —^.  (7) 

r 

J,  is  the  first-order  Bessel  function.    The  arguments  of  the  sombrero  functions  are  the  distances 

^1  =  ^ix-wf+(y-z)^. 


r,  =  ^{w-af^(J-b)\ 


where  r,  is  the  distance  between  two  points  in  the  object  plane,  rj  is  the  distance  between  the  first  point 
in  the  object  plane  and  the  point  in  the  image  plane,  and  rj  is  the  distance  between  the  second  point  in 
the  object  plane  and  the  point  in  the  image  plane. 

The  first  sombrero  fiinction  is  the  result  of  the  of  the  Van  Cittert-Zemike  theorem.  The  second 
and  third  sombrero  functions  are  both  the  Green's  function  of  the  microscope.  We  normalize  Equation 
(6)  by  dividing  by 

N(fl,b)  =  C\dx\dy\dw\dzsomb{air^somb(r^somb{r^,  /m 

a  -  ^  -  0, 

which  is  the  intensity  at  the  origin  in  the  absence  of  the  half-plane.    The  result, 

/      =    ija^b)  (10) 

"""       /V(0,0)' 


nor 


is  a  function  that  asymptotically  approaches  0  as  y  approaches  -oo  and  approaches  1  as  y  approaches 

+00. 
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The  resulting  integrals  cannot  be  solved  analytically.  We  evaluated  these  integrals  numerically 
for  several  coherence  parameters  and  computed  the  edge  profiles,  intensity  as  a  fiinction  of  distance 
along  the  y  axis,  for  several  coherence  parameters  (see  Appendix). 

Figure  1 5  shows  edge  profiles  calculated  for  a  number  of  coherence  parameters. 
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Figure  15.         Edge  profiles 

1  resolution  limit  =  3.14. 

We  also  compiled  a  table  of  the  edge-finding  intensities  for  several  coherence  parameters. 
Figure  16  is  a  graph  of  those  data.  We  can  interpolate  for  any  value  of  the  coherence  parameter 
between  0  and  2.  The  edge  intensity  starts  at  0.25  for  coherent  light,  a  =  0,  and  approaches  0.5  for 
incoherent  light,  a  =  oo. 
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Figure  16.         Edge  intensity  as  a  function  of  coherence  parameter. 

Because  the  edge-finding  algorithm  in  the  computer  program  uses  a  relative  intensity,  we  need 
to  know  the  numerical  aperture  of  the  illumination  spot.  We  know  the  dimension  of  the  spot  of  light, 
but  we  need  to  measure  the  distance  from  the  end  of  the  fiber  to  the  metal  plate.  We  did  this  by  using 
the  micrometer  for  coarse  focus.  We  retracted  the  fiber  into  the  SRM  box,  moved  the  white  plate  until 
it  touched  the  cover  slip,  and  noted  the  position  on  the  micrometer.  We  then  backed  off,  moved  the 
fiber  out  of  the  SRM  box,  focused  on  it,  and  noted  the  position.  The  difference  between  the  two 
positions  is  the  distance  between  the  spot  and  the  objective.  We  need  to  find  the  distance  to  the  focal 
point  of  the  objective,  however,  not  to  the  surface  of  the  lens.  Since  the  working  distance  (the  distance 
between  the  front  surface  of  the  lens  and  the  focal  point)  of  the  objective  is  0.6  mm,  we  subtracted 
0.6  mm  to  find  the  distance  between  the  spot  and  the  fiber  end,  or  5.47  mm.  When  we  moved  the  plate 
to  the  cover  slip,  it  met  a  gradual  resistance,  not  a  sudden  stop.  We  estimate  that  the  point  at  which 
the  plate  just  touches  the  cover  slip  is  within  0.25  mm;  this  is  the  estimated  3a  uncertainty. 

Using  the  3a  uncertainty  of  the  distance  between  the  spot  and  the  fiber  end.  we  estimated  the 
3a  uncertainty  in  the  mean  diameter,  noncircularity,  and  major  axis  angle.  We  then  divided  by  3  to 
estimate  the  la  uncertainty  of  the  coherence  parameter:  the  coherence  parameter  is  1.27  with  a  3a 
uncertainty  of  ±0.02  because  of  the  uncertainty  in  the  fiber-to-spot  distance.  Using  the  value  1.27,  we 
estimate  the  theoretical  edge-finding  intensity  as  0.414.    We  estimated  the  uncertainty  of  the  edge- 
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finding  intensity  by  calculating  the  edge  intensity  for  coherence  parameters  of  1.25  and  1.29;  these 
differ  from  1.27  by  ±0.02.  We  thus  estimate  the  3a  uncertainty  of  the  edge-finding  intensity  as  ±0.02. 
We  did  this  calculation  before  taking  any  data  for  this  paper,  and  we  used  the  0.414  intensity 
for  finding  all  edges.  The  3ct  uncertainty  in  the  mean  diameter  corresponding  to  the  3cj  uncertainty  in 
the  edge  intensity  is  5  nm  and  is  insignificant  compared  to  the  uncertainty  due  to  defocusing.  The  la 
uncertainty  in  the  edge  is  1.7  nm.  The  la  uncertainty  in  the  noncircularity  is  ±0.000  07%;  the  la 
uncertainty  in  the  major  axis  angle  is  ±0.07°. 


7.    Simulations 

We  used  an  edge  profile  from  Figure  16  to  run  a  simulation  of  how  the  pixels  in  the  camera 
affect  the  image.  While  we  were  looking  at  plots  of  the  polar  edge  tables,  we  noticed  that  certain 
aspects  of  the  edge  table  did  not  shift  in  angle  as  the  fiber  was  rotated.   Figure  1 7  is  a  polar  edge  table. 
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Figure  17.         Edge  table  of  a  highly  circular  fiber. 

The  edge  table  data  display  an  approximately  sinusoidal  behavior  as  well  as  what  appears  to  be 
random  noise.  When  the  fiber  is  rotated,  the  maxima  and  minima  of  the  distribution  change  position; 
this  indicates  that  they  result  from  the  noncircularity  or  cleave  angle  of  the  fiber. 
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Figure  18  shows  the  same  data  as  in  Figure  17  but  expanded  around  1.57  rad  (90°).  On  this 
scale,  the  data  display  nonrandom  oscillations  at  1.57  rad  and  4.71  rad  (270°).  These  nonrandom 
oscillations  do  not  change  position  when  the  fiber  is  rotated,  so  they  are  not  related  to  the  geometry  of 
the  fiber. 
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Figure  18.         Oscillation  of  the  image  of  the  fiber  edge  . 

Since  these  nonrandom  areas  do  not  move  with  the  fiber,  they  must  be  some  artifact  of  the 
microscope.  We  used  the  theoretical  edge  profile  data  for  a  coherence  parameter  of  1.28  to  simulate 
the  effect  of  finite,  rectangular  pixels  on  the  image  of  a  circular  fiber.  We  approximated  the  continuous 
two-dimensional  edge  profile  of  the  fiber  as  a  very  fine  grid  of  discrete  points.  The  intensitv'  of  each 
point  was  calculated  from  an  edge  profile  based  on  the  distance  of  the  point  from  the  center  of  the  fiber. 
We  approximated  the  pixels  as  rectangular  regions  of  uniform  responsivity.  The  intensity  recorded  b> 
one  of  the  simulated  pixels  is  the  sum  of  the  values  at  the  discrete  points  divided  by  the  number  of 
points  within  the  pixel.  Since  the  fiber  edge  curves  only  negligibly  over  a  distance  equal  to  the 
resolution  limit  of  the  objective  lens,  we  can  use  the  edge  profile  for  a  half-plane  to  approximate  the 
edge  profile  of  the  curved  fiber.  In  this  approximation,  the  intensity  as  a  function  of  the  fiber  radius 
is  the  same  as  the  intensity  along  the  y  axis  for  the  half-plane  calculation.  We  assumed  a  perfectlN 
circular  fiber.    Figure  19  is  the  edge  table  of  the  simulated  fiber  image. 
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Figure  19.         Simulated  edge  of  a  circular  fiber. 

Figure  19  shows  that  a  circular  fiber  imaged  with  rectangular  pixels  theoretically  appears 
noncircular  with  a  noncircularity  of  0.025  %  ;  the  major  axis  is  horizontal.  This  apparent  noncircularity 
is  significant  compared  the  0.05  %  noncircularity  of  the  fiber  we  have  been  using. 

The  radius  we  chose  for  the  simulation  was  62.5  \xm.  Figure  19  shows  a  0.015  |im  offset 
between  the  apparent  mean  radius  and  the  true  radius.  While  this  is  significant,  it  is  a  part  of  the  offset 
between  the  micrometer  measurement  and  the  average  of  the  defocus  runs  and  does  not  represent  a 
source  of  error. 

Figure  20  is  the  simulated  edge  table  expanded  about  1.57  rad  (90°).  Figure  20  has  nonrandom 
oscillations  at  1.57  rad  (90°)  and  4.71  rad  (270°)  similar  to  those  in  the  real  data  (Figure  18).  The 
magnitude  of  the  nonrandom  oscillations  predicted  by  the  simulation  is  smaller  than  what  we  actually 
observed  in  the  real  data.  Further,  the  simulation  predicts  a  similar  result  at  0  rad  and  3.14  rad  (180°). 
We  do  not  see  this  in  the  actual  data  at  those  angles. 
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Figure  20.         Oscillations  of  the  image  of  a  simulated  edge. 


The  simulation  also  predicts  an  overall  sine  curve  with  maxima  at  0°  and  180°  (Figure  19)  and 
a  noncircularity  of  0.025  %.  We  do  not  see  this  in  the  real  data  either.  Since  yaw  produces  an  effect 
just  opposite  to  what  is  predicted  by  the  simulation,  a  certain  yaw  angle  could  exactly  cancel  the  overall 
sine  due  to  pixelation.  Since  we  adjusted  the  tilt  to  remove  any  sine  curve  which  did  not  rotate  with 
the  fiber,  we  may  have  inadvertently  corrected  for  the  sine  curve  due  to  the  pixelation. 

The  difference  in  magnitude  between  the  simulated  and  real  nonrandom  oscillations  indicates 
that  we  are  not  taking  something  into  account.  We  looked  at  the  data  from  the  camera  to  check  the 
shape  of  the  intensity  profile.  Figure  21  compares  the  theoretical  and  real  edge  profiles  (the  position 
of  the  real  profile  on  the  horizontal  scale  is  arbitrary).  The  edge  profile  of  the  real  data  is  broader  than 
that  of  the  theoretical  calculation. 
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Figure  21.         Theoretical  and  real  edge  profiles 
1  resolution  limit  =  3.14. 

The  difference  between  the  shapes  of  the  measured  and  theoretical  edge  profiles  probably 
accounts  for  the  difference  in  the  magnitudes  of  the  nonrandom  oscillations.  It  probably  also  means  that 
the  intensity  we  have  been  using  for  edge  detection  is  incorrect;  this  could  cause  all  or  part  of  the 
0.19  |j.m  error  in  the  mean  diameter. 

The  measured  edge  profile  is  probably  broader  than  the  theoretical  because  of  aberrations  in  the 
objective  lens.  Our  calculation  did  not  include  lens  aberrations,  which  increase  the  radius  of  the 
Green's  function  of  the  system  and  therefore  change  the  mutual  coherence  function.  In  addition,  we 
assumed  a  mutually  incoherent  source,  even  though  the  white  plate  was  illuminated  by  a  source  that  was 
itself  not  wholly  incoherent.  These  factors  make  our  calculation  quantitatively  incorrect  and  defeated 
our  attempt  to  measure  the  diameter  of  the  fiber  without  recourse  to  the  calibration  artifact. 

One  of  the  final  stages  in  the  camera  is  a  low-pass  filter.  This  low-pass  filter  operates  on  the 
electronic  signal  before  it  gets  to  the  frame  digitizer.  This  filter  sees  the  outgoing  data  as  a  long  string 
of  voltage  peaks.  The  net  effect  is  a  low-pass  filtering  of  each  horizontal  line  of  pixels  but  not  of  the 
vertical  columns  of  pixels.  We  assumed  that  the  low-pass  filter  is  a  simple  resistor-capacitor  filter 
whose  step  response  is  an  exponential  decay  with  time.    The  cutoff  frequency  of  the  filter  corresponds 
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Figure  30.         Mean  diameter  as  a  function  of  fiber  rotation  angle. 


Figure  31  shows  noncircularity  as  a  fiinction  of  fiber  rotation  angle.  The  noncircularity 
oscillates  about  the  micrometer  value  as  the  fiber  rotation  angle  changes.  This  could  indicate  that  the 
beam  splitter  mount  has  moved  and  created  a  tilt  angle. 
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Figure  31.  Noncircularity  as  a  function  of  fiber  rotation  angle. 

Figure  32  shows  the  angle  of  the  major  axis  as  a  function  of  fiber  rotation  angle.     The 
agreement  of  these  data  is  very  good. 
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Figure  32.         Angle  of  major  axis  as  a  function  of  fiber  rotation  angle. 
12.    Conclusions 

With  an  SRM  fiber  as  a  calibration  artifact,  we  succeeded  in  making  mean  diameter 
measurements  with  an  uncertainty  less  than  0.1  |j.m.  However,  we  did  not  succeed  in  making  the 
measurements  absolutely,  that  is,  without  using  the  calibration  artifact. 

The  difference  between  a  real  edge  profile  and  the  theoretical  edge  profile  in  Figure  22  shows 
that  we  have  not  calculated  the  edge-finding  intensity  in  a  realistic  manner.  This  is  probably  the  result 
of  aberrations  in  the  lens,  which  was  an  inexpensive  microscope  objective.  A  better  lens  might  reduce 
the  error  to  less  than  the  uncertainty  of  the  measurements.  Alternatively,  the  calculation  of  the  edge 
profile  could  be  altered  to  take  aberrations  into  account. 

The  3a  uncertainty  of  0.07  )im  in  our  measurements  of  mean  diameter  could  be  reduced  by  a 
better  focusing  method.  Ideally,  we  would  like  computer-controlled  motorized  micrometers  to  move 
the  fiber.  The  computer  could  continuously  adjust  the  micrometers  to  produce  the  sharpest  edge 
possible.  Motorized  micrometers  alone  might  improve  focusing  simply  because  placing  our  hands  on 
the  instrument  moves  the  fiber  noticeably.  Likewise,  the  3a  uncertainty  of  0.03%  in  the  noncircularity 
could  be  improved  by  motor-driven  tilt  adjustments  with  computer  feedback. 

Finally,  details  of  the  software  used  to  create  the  edge  tables,  fit  the  ellipses,  calculate  the 
centers,  and  so  on  are  provided  in  [13]. 
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Appendix. 
Evaluation  of  the  Integrals  of  Partial  Coherence 


Problem  Statement 


A  uniform  disk  source  of  incoherent,  quasimonochromatic  light  obstructed  by  a  half-plane  and 
focused  by  a  lens  system  with  circular  stop  produces  image  plane  intensity  at  a  point  (a,  b)  modeled  using 
the  van  Cittert-Zemike  approximation  [9],  [10]  as 

/cxD  /-oo  /"oo  roo 

dx         dy  dw         dz  s{ari)  s{r2)  sirs)  e'^^'''     '''\  (1) 

-oo         Jo  J-oo         Jo 

where  the  sombrero  function  s  is  defined  by  the  formula 


s{r)  =  -^^,  (2) 


Ji  is  the  Bessel  function  of  order  1 ,  and 


n  =  ^{x-w)^  +  {y-zy,  (3) 


r2^^J{x-ay  +  {y-b)^  (4) 

ra  =  yjiw  -  a)2  +  [z  -  6)2,  (5) 

r4  =  ^x^+y^,  (6) 

rs  =  yw'^  +  z'^.  (7) 


The  intensity  at  (a,  6)  without  the  half-plane  obstruction  is  similarly  modeled  as 

/oo  /-oo  roo  fOO 

dx  dy  dw  dz  s{ari)  s{r2)  sirs)  e'^^^'  -^'  \ 

-oo  J —oo  J  —  OO  V— 00 

and  the  relative  intensity  is  /(a,  b)/Nia,  b). 


(8) 


The  constants  a  and  /3  are  determined  by  the  geometry  of  the  experimental  apparatus.  For  a  point 
source  of  (coherent)  light,  a  —  P  =  0,  and  the  integrals  separate.  In  the  experimental  setup  of  interest,  a 
had  a  value  near  1  and  /5  was  roughly  10""*.  The  relative  intensity  was  desired  for  a  =  0  and  a  range  of  real 
values  for  b. 

The  computation  of /(a,  b)  using  standard  software  would  be  prohibitively  expensive.  We  describe 
several  simplifications  which  reduce  the  number  of  operations  enough  to  make  the  computation  possible  on 
a  workstation. 

Numerical  Method 

We  observe  that  for  small  A6,  /(a,  b)  and  /(a,  b+Ab)  differ  by  an  integral  that  is  easier  to  compute 
than  either  integral  for  /  alone.  More  precisely,  9/(a,  b)/db  is  an  integral  in  three  dimensions  (rather  than 
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four): 


—/(a,  6)   =   \im6-^{I{a,b  +  S)-I{a,b)) 
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The  same  method  can  be  used  to  show 


—Nia,b)  =  -Nia,b)=0, 


(9) 


(10) 


so  A'^  depends  only  on  a  and  /5. 

Simpson's  quadrature  rule, 

-b+h 


ro+n  u 

/         fix)  dx^-  ifib)  +  4/(6  +  h/2)  +  fib  +  h))  , 
Jb  6 


(11) 


is  employed  to  compute  an  approximation  /(a,  6)  to  the  intensity  /(a,  6): 

7(a,  6  -f  /i)  =  /(a,  6)  +  ^  (^/(a,  b)  +  4^/(a,  6  +  h/2)  +  ^/(a,  6  +  /i)  )  , 
for  a  chosen  step  size  /i.  This  approximation's  error  is  of  order  Oih'^)  (fourth  order). 

We  approximate  the  integral  (9)  for  dl/db  by  restricting  the  domain  to  be  finite, 

rt  rt  rt 


(12) 


It  =  2        dx       dy        dw  siari)  sir2)  sirs) 
J-t      Jo       J-t 


APir^'-rs') 


(13) 


2  =  0 
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which  allows  the  integration  to  be  performed  to  high  accuracy.  To  obtain  an  estimate  of  /qo,  we  make  the 
extrapolation  assumption  that  It  has  the  form 

We  truncate  the  right  side  of  (14)  at  three  terms,  obtain  3  linear  equations  in  the  3  unknowns  /qc,  «!,  and 
a2  by  substituting  computed  values  for  It,  ht,  and  l4t,  and  then  solve  for  the  unknowns.  This  extrapolation 
procedure  can  be  expected  to  be  reliable  if  t  is  chosen  large  enough  that  It  —  ht  is  small  relative  to  It  and 
It  —  ht  ~  2(/2f  -  ht)-  Our  computations  indicate  that  these  conditions  are  met  for  the  parameter  values 
of  interest,  a  ^  1,  P  ^  10""*,  a  =  0,  and  — 6  <  6  <  6,  provided  that  t  >  20.  Based  on  the  computations, 
we  believe  that  this  extrapolation  decreases  the  error  in  approximating  /qo  from  0{t~^)  to  0(i"~'^);  that  is,  it 
yields  third-order  convergence. 

We  compute  h  with  Gaussian  quadrature  [14],  [15].  In  particular,  we  use  Gauss-Legendre  quadra- 
ture with  n{t)  nodes  in  y  and  2  n{t)  nodes  in  x  and  w,  where 

n(t)=\\      "11^''  (.5) 

I    at,     a  >  I. 

Gaussian-Legendre  quadrature  approximates  an  integral  by  an  linear  combination  of  values  of  the  integrand, 

■1 


/_ 


f{x)dx  ^^Wrf{Xi),  (16) 


1  ^=1 


where  the  nodes  xi, . . .  ,Xn  and  weights  iwi , . . . ,  lo^  are  determined  so  that  (16)  is  exact  if  /  is  a  polynomial 
of  degree  up  to  2n  —  1.  (The  nodes  coincide  with  the  roots  of  the  Legendre  polynomial  of  degree  n.)  We 
determine  the  nodes  and  weights  using  standard  software  [16].  These  quadratures  converge  extremely  rapidly 
for  any  infinitely  continuously  differentiable  function  /;  in  our  case  the  quadratures  are  accurate  to  machine 
precision  (15  digits). 

The  integrand  of  It  involves  the  Bessel  function  Ji,  which  we  compute  from  [17] 

2n 

Jn-l{x)  =  Jn{x)  -  Jn+l{x),  (17) 

X 

applied  in  order  of  decreasing  n,  combined  with  the  normalization  equation 

1  =  Jo(x)  +  2h{x)  +  2J4{x)  +  2Mx)  +  ■■■.  (18) 

Since  the  values  of  s{x)  =  Ji{x)/x  are  required  repeatedly,  we  tabulate  this  function  at  equispaced  values 
of  X-  from  0  to  the  maximum  value  required  (\/5  ^  max{a,  1})  and  obtain  the  function  at  required  values  of 
X  by  fourth-order  interpolation.  The  table  is  set  up  so  that  the  interpolation  operation  requires  only  three 
multiplications  and  three  additions.  The  tabular  increment  of  x  of  0.01  is  chosen  to  yield  10  digit  accuracy. 

A  straightforward  computation  of  I  {a,  b)  proceeding  directly  from  its  definition  (1)  would  be  be- 
yond our  current  computational  resources.  The  result  of  our  optimizations  is  that  the  required  integrals  can 
be  computed  on  a  workstation.  We  use  a  high-performance  workstation  that  has  peak  performance  of  over 
50  million  floating-point  operations  per  second;  for  a  given  value  of  a,  /5,  and  a,  we  obtain  I{a.  h)  for  -6  < 
6  <  6  in  under  an  hour  of  computation  time.  The  resulting  accuracy  is  limited  by  the  accuracy  of  the  extrap- 
olation; for  t  =  20  the  relative  error  is  less  than  10""^. 
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Supplementary  Analytical  Result 

The  computed  results  for  /3  ^  10~^  nearly  coincide  with  those  for  /3  =  0;  the  fact  that  this  holds 
for  the  integrals  can  be  shown  analytically.  As  asserted  above,  N{a,  b)  is  actually  independent  of  (a,  b),  and 

for  (3  =  0,  N{a)  =  N{a,  b)/{27T^)  has  the  values 


1,         0  <a  <  1 
a~^,     Q  >  1, 


N{a)  =  {    l'_2      :""-'  (19) 


which  Ron  Wittmann  has  demonstrated  analytically  [18]. 
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